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Abstract 

The nonlinear evolution of long wavelength non-stationary cross-flow vortices in a 
compressible boundary layer is investigated and the work extends that of Gajjar (1994) to 
flows involving multiple critical layers. The basic flow profile considered in this paper is 
that appropriate for a fully three-dimensional boundary layer with 0(1) Mach number and 
with wall heating or cooling. The governing equations for the evolution of the cross-flow 
vortex are obtained and some special cases are discussed. One special case includes linear 
theory where exact analytic expressions for the growth rate of the vortices are obtained. 
Another special case is a generalization of the Bassom & Gajjar (1988) results for neutral 
waves to compressible flows. The viscous correction to the growth rate is derived and it is 
shown how the unsteady nonlinear critical layer structure merges with that for a Haberman 
type of viscous critical layer. 



1. Introduction 


The transition to turbulence of fluid flows has aroused considerable interest in recent 
years and substantial progress has now been made towards understanding the early stages 
of transition in many shear and planar boundary layer flows. The use of self-consistent 
asymptotic analysis in conjunction with the correct description of the flow structure in 
the limit of large Reynolds numbers, has given a remarkable insight into the delicate 
balances and different processes operating in a transitional boundary layer. There are 
many routes to transition in a boundary layer. Carefully controlled experiments with pre- 
dominantly two-dimensional disturbances have highlighted the role played by disturbance 
three-dimensionality, Nishioka et al (1979), Klebanoff et al (1962), Kachanov Sz Levchenko 
(1984), Corke & Mangano (1989). After an initial stage in which the two-dimensional 
disturbances dominate, the onset of three-dimensionality leads rapidly to more complex 
structures and nonlinear interactions leading ultimately to turbulence. Much of this initial 
linear phase and some of the subsequent complicated nonlinear stage can now be described, 
if not quantitatively at least in a convincing qualitative manner using asymptotic theo- 
ries, see for example the papers by Goldstein & Choi (1989,1990), Hall & Smith (1991), 
Goldstein & Lee (1992), Wu et al (1993,1994), Smith, & Stewart (1987), Stewart & Smith 
(1992), Smith & Bowles (1992), Hultgren (1992), Goldstein (1994). 

Our concern here is with the stability of fully three-dimensional boundary layers rather 
than two-dimensional boundary layer flows. One of the main differences between the 
stability of two and three- dimens ional incompressible boundary layers is the occurrence 
of cross-flow instability in the latter. Whereas a typical two-dimensional incompressible 
boundary layer is not subject to invisdd instability, a three-dimensional boundary layer 
almost invariably is. The mechanism was first explained by Stuart in Gregory, Stuart 
&: Walker (1955) and he showed that in certain flow directions the velocity profile could 
be inflexional and hence subject to Rayleigh instability. This instability is known as 
cross-flow instability. A two-dimensional compressible boundary layer is also subject to 
inviscid instability because of the generalized inflexion point, Lees & Lin (1946), Mack 
(1984, 1986). However cross-flow instability is believed to be more important even for 
compressible flows, especially at low Mach numbers, because of the much larger growth 
rates. Cross-flow instability is the main cause of transition of several flows of aerodynamic 
importance such as the flow past swept aircraft wings and the flow near the attachment 
line. See for instance the reviews of three-dimensional boundary layer instability by Reed 
& Saric (1989) and Arnal (1986). Despite its importance cross-flow instability has not 
received as much theoretical attention as the planar case as it deserves. 

Many experiments have been performed to study cross-flow instability. The flow over 
a rotating disk has been widely used as a prototype of a fully three-dimensional bound- 
ary layer flow and has been studied experimentally by Gregory, Stuart & Walker (1955), 
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Wilkinson & Malik (1983) amongst others. The most commonly observed manifestation 
of cross-flow instability is the formation of a pattern of equi-spaced co-rotating stationary 
vortices. The stationary cross-flow instability arises from the amplification of zero fre- 
quency waves in the boundary layer. A non-stationaxy, or a travelling wave, pattern of 
cross-flow vortices is also possible. In experiments on swept- wing flows Bippes& Nitchke- 
Kowsky(1987), Muller & Bippes (1988), Kohama et al (1991), both types of instability 
have been found to occur. Poll (1985) has also observed both types of instability in ex- 
periments with swept cylinders. Muller Bippes (1988) have also investigated nonlinear 
effects and found that both the stationary and the non-stationaxy vortices reach a non- 
linear eq uili brium amplitude. In the experiments of Kohama et al (1991) an explosive 
secondary instability was seen prior to the onset of transition to turbulence. 

Theoretical investigations of the linear stability of cross-flow vortices have to date 
been largely confined to incompressible flows. The stability properties of stationary cross- 
flow vortices were first investigated by Gregory, Stuart & Walker (1955). Malik(1986) and 
Balakumar &: Malik(1990) have computed neutral and growth rate curves for both the 
stationary and travelling wave disturbances for the rotating disk. Hall (1986) computed 
the corrections to the neutral values for the stationary case in the in viscid limi t and also 
obtained the structure for the zero shear stress mode discovered by Malik (1986). Bassom 
& Gajjar (1988) have looked at various linear and nonlinear aspects of the neutral stability 
of non-stationary cross-flow vortices. There have been various studies based on numeri- 
cal simulation of either the full Navier- Stokes equations, Meyer & Kleiser (1988) or an 
approximate form such as parabolised stability equations approach of Malik & Li(1992), 
Malik it et al. (1994). The numerical simulations show the nonlinear saturation of cross- 
flow vortices which is in qualitative agreement with the experiments of Muller & Bippes 
(1988). 

An objective of the present work is to study the nonlinear evolution of non-stationary 
cross-flow vortices in compressible boundary layers using unsteady no nlin ear critical layer 
theory. The latter has been successfully used in the studies of the instability of m a n y planar 
boundary layer and shear flows, see for instance the papers cited earlier. Excellent reviews 
of the early ideas on nonlinear critical layers are given in Stewart son (1981) and Maslowe 
(1986), and more recent work is surveyed in Goldstein (1994) and Cowley Sz Wu (1994). 
The assumptions and scalings used here lead to a fully nonlinear partial differential system 
for the evolution of the vortex as in Goldstein, Durbin & Leib (1987), Goldstein & Leib 
(1988), Goldstein & Hultgren (1988), Goldstein & Wundrow (1990), Gajjar (1994), and 
not an integro- differential equation of the Hickemell (1984) type. The scalings adopted 
are appropriate for low frequency long wavelength cross-flow vortices which have small 
growth rates. For shorter waves the weakly nonlinear approach leading to a Hickernell 
type amplitude equation is more appropriate, Gajjar (1995). Gajjar (1994) has studied the 
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development of long wavelength oblique modes on two-dimensional compressible boundary 
layers and the present work extends this to three-dimensional boundary layers. The bulk of 
the analysis can be simplified through the introduction of appropriate Squire coordinates 
and provided variations in the ‘cross-Squire’ directions are small the analysis is similar to 
the two-dimensional case. The situation where the cross-Squire variations axe not small 
has been examined in Gajjar & Arebi (1995) for non-stationaxy vortices, and in Gajjar 
(1995) for the stationary cross-flow vortex. In both cases the analysis lead to a novel 
integro- differential equation for the vortex amplitude involving spanwise derivatives of the 
amplitude in the kernal function as in Wu (1993). 

For low frequency waves, for certain wave angles, there exists two critical layers in the 
flow one located near the wall and the other in the main part of the boundary layer. The 
flow structure and analysis axe therefore very similar to Gajjar (1994) with modifications 
for the additional critical layer which occurs in the three-dimensional case. Implicit in the 
present analysis is that there is no region of locally supersonic flow near the wall, in the 
terminology of Mack (1984). This places a certain restriction on the Mach number. 

In section 2 below the problem is formulated for a compressible fluid satisfying the 
perfect gas law, and in sections 3 through 5 the analysis leading to the equations governing 
the evolution of the vortex is given. Finally we end with some discussion in section 6. 

2 Formulation of the Instability Equations 

Consider coordinates x* = (* t ,y t ,« t ) such that x t is directed along the wall, y t is 
perpendicular to the wall and z 1 is in the spanwise direction and such that the corre- 
sponding velocity components axe given by ut = (ut, vt,ti>t) } where a suffix f denotes a 
dimensional quantity. The Navier- Stokes equations axe non-dimensionalised by taking L 
to be a typical streamwise lengthscale and 17^,, a free stream velocity. Also a quantity 
(frlo below denotes a local free stream value. The coordinates (x\y\z^) axe written as 
( Lx,Ly,Lz ), the velocity (u t ,v t ,w t ) as (U^u, U^v, U^w), the pressure is scaled on 
the viscosity pt is scaled on /x^p, the bulk viscosity At is scaled with respect 
to pt^A, the temperatute T* is scaled with respect to T^T, the density pt is scaled with 
pt^p, and finally the conductivity fct is scaled with k^k. The non-dimensional form of the 
equations axe: 

t = °’ (2 ' la) 


dui du{ 1 d 

dt dxj p dxj 


(P*ij) + + ( A “ §/*)«**)}» 


R I dxj 


dp dp 
— 4- Ui — — 
dt l dxi 


= T P 


+ 


dp dp 

J_ Xfi 

dt ^ X dxi 
o d 


(r-i) fi 

r 


R 


— — J. / \ 2 i / OH % 

+ p- o ~ I o / ie *j c v + (^ — 


. du{ . 
dxi 


RM^dx 


i \PTdXi) 


(2.16) 


(2.1c) 
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TMlp = P T, 


(2.1 d) 


( dui duj \ 

6,-7 l dxj dx{ / 

In the above equations the Reynolds number R, Mach number Moo, 
number Pr axe defined by 


(2.1e) 

and the Prandtl 


R = 


ULLpl 


l^oo 


M 2 = 

CO 


u£ 

TUTo 


Pr = ^ 

kl 


t 


and c p is the specific heat at constant pressure, T is the ratio of specific heats, and R is 
the gas constant. The Reynolds n um ber is assumed to be large and the Prandtl number 
to be constant. 

The leading order three-dimensional basic boundary layer flow is taken to be of the 
form 


U s = (U b (x,Y,z),R-^ 2 V b (x,Y,z),Wb(x,Y,z)), (2.2a) 

P=Pb(*, z ), T = Tb(x,Y,z), PB = Pb( x ,Y,z), (2-2 b) 

where y = R~*Y defines the boundary layer coordinate Y . 

The analysis is simplified through the introduction of Squire coordinates (£, z), see 
Figure 1, along and normal to the direction of propagation of the wave, as follows 

£^ RHah ^ +a e±-^), (2.3a) 

±^Ri(/3 h ~ + l3h^ + a h ~), (2.3b) 

I- - Rih.*(-ac + (2-3c) 

dt oQ 


where the quantities z and X in (2.3) are defined by £ — R*h(ax + /? z achi), 
z = R*h(az - (3x) and X = R*h 2 (ax + (3z). The wave is taken to grow on the slow 
X scale. In equations (2.3) a and (3 are defined by a = cos#, = sin#, where # is the 
direction of propagation of the wave. Furthermore c denotes the phase speed of the wave 
and <t is a real quantity. The quantity h is a scale factor (h << 1) introduced such that 
we are considering long waves (based upon the Rayleigh scalings) and the growth rates are 
therefore 0(h 2 ). 

The cross-flow velocity field is defined by u = au + f3w, w = — flu + aw with respect to 
the chosen Squire coordinates. With respect to these coordinates the basic flow velocities 
are given by, 

U B = cJJ B +pW B , W B = -flU B + aW B , 
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in the £ and z directions respectively. 

The stability analysis depends crucially on the properties of the basic flow near the 
wall and the upper critical level and it is assumed that 

U B ~ X x Y + A 2 F 2 + ... as Y^0+, 

PB ~ Ro 4- RiY + . . . as Y — ► 0+, 

Tb ~ Sq -f- •S’iT -f- • • • as Y — * 0+, 

W b ~6 1 Y + 6 2 Y 2 + ... as Y — + 0+, 

Ub ~ X^Y -Yu) + \' 2 (Y -Yu) 2 + ... as Y — > Yu, (2 ‘ 4) 

PB ~ JRJ + - Yu) + . . . as Y ^ Yu, 

T b ~ 5* + S*(y -Yu) + ... as y - Y v , 

W b ~S;+6;(Y -Yu) + ... as T-Ty, 

where the upper critical level is defined by Ub — 0 at Y = Yu. In addition at the outer 
edge of the boundary layer we require that 

U B ~ Uoo, W B ~ 1^005 PB ~ Poo, T b ~ Too, as Y — * oo, 

and these outer edge values may depend on x and z in general. 

A typical cross-flow velocity profile is of the form as shown in Figure 1 so that Ai and 
AJ are of opposite signs. Furthermore the combination Ajc is positive. 

The main purpose of the analysis below is to determine the wavenumber and the linear 
and nonlinear growth rate of the vortex for a presribed mean flow and given frequency u>. 
The wave-angle 0 is taken to be a parameter and for a fixed value of 6 the position Yu is 
also fixed. 

3 Stability analysis outside critical layers. 

We carry out a perturbation analysis by introducing small disturbances of size 8 and 
considering expansions for the perturbed flow quantities. If we ignore the upper critical 
layer the bulk of the expansions below follow from those given in Gajjar (1994). The 
disturbance size 6 will be fixed subsequently to be of 0(/i 3 ) such that nonlinear effects first 
alter the linear growth rate. For clarity however, and since the flow structure outside the 
critical layer is largely dictated by linear dynamics, we work with 8 for the time being. We 
first consider the main part of the boundary layer flow for which Y = 0(1) and expand 
the various flow quantities as 

u = U B + 8{u i + hu 2 + •••)> 

v = 8(hvi d- h, 2 v 2 -J- . . .), 
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(3.1a) 

(3.16) 



W = Wb + 6(wi + hw2 + ...), 

(3.1c) 

P = PB + 6(hpi + h?p2 + ...), 

(3.1(f) 

p = PB + &{pi + hp 2 + . . •), 

(3.1c) 


(and the perturbed quantities are taken to be independent of z to this order). 

These expansions (3.1) are substituted into the governing equations (2.1) yielding the 
following solutions at leading order 


Hi = AoUb 


Y 5 


Vl 


= - Ao ( Ub , 


(3.2a, b) 


pi = Pi, pi = AopByi — Ao(Wb y (3.2c, d, e) 

where Aq is given by Ao = A(X)exp (17^) + c - c * -Pi = Ai(-^0 ex P (*7C) + c - c - 7 is 
the wavenumber of the cross-flow vortex. 

At next order the analysis yields 


P 2 = P 2 ~ 



UbPBVi ( <IY, 


(3.3 a) 


V2 


=u B r — 

Jy 0 PB 


Pl hrdY + acA 0< 


-AfU B - 


PbU 2 b 


UbPi< v 

T Pb ‘ 


W 2 


27 




(3.36) 


(3.3c) 


The constant Yq is chosen to be such that Yo > Yu* It can be seen that the integrals in 
equation (3.36c) are singular at the locations where Us — 0 and this occurs at Y = 0 and 
Y = Yb (see Figure 2). Additional expansions must be introduced in order to obtain valid 
expressions for the flow quantities near Y = 0 and Y = IV* I* 1 the neighbourhood of the 
lower critical layer situated near the wall we introduce regions Jrfc of thickness 0{K) and in 
the vicinity of Yu (the upper critical level situated in the m ain part of the boundary layer 
flow) regions IVii are introduced (again of thickness O(h)). Hence the Y = 0(1) region 
is divided into two regions /Ji due to the presence of the upper critical layer Yu in the 
main part of the boundary layer, see Figure 2. Moreover we define the unknown functions 
in equation (3.36) such that A ^ relates to region JJ— and A ^ relates to region //+ 
where A± (and P 2 ) are functions of £ and X. 


Considering first the regions I± of thickness 0(h) close to the wall we introduce a 
coordinate Y such that Y = kY and the above solutions therefore imply 


u == h\\Y -h h 2 \ 2 Y 2 + . * • "I - 6(ui “h hv, 2 "I - • • •)> (3.4u) 


7 



v — 6(h 2 vi + h z v 2 + ...), (3.4 6) 

P = PB + f>(hpi + h 2 p 2 + ), (3.4c) 

P — Ro + hR\Y + h 2 R 2 Y 2 + S(pi + hp 2 + ...), (3.4<i) 

w = So 4- hS\Y + . . . + + hw 2 + ...). (3.4e) 


Substituting into the governing equations we obtain at leading order the solutions 


ui = = 

pi = Pi exp(t7C), pi ( = 


~P i< 


- A 0( (A!y - ac), 
—viRi _ 8\vi 


(A xY -ac) 


wi( = 


(\iY-a c) 


(3.5 a, b) 
(3.5 c,<£) 


We note the density function has a pole singularity at lx, = ac/Xi (as indeed does the w\ 
fluctuation). The requirement that Vi — ► 0 as Y — ► 0 gives 


Pi — A olcXiRq. 


(3.6) 


At the next order we find 


v 2 = (—aA 0( — acA 0x ) — 


X 2 A 0( (t} 2 - a 2 c 2 ) 

A? 


+ ^aT^‘ RiO£2c2 ^ 0c _ ^ 2< 


acA 0( 

Ai 


+ W + tth, 

P2 =p2((,X), 


(3.7 a) 
(3.76) 


where <f > ^ are terms arising from the continuation across the lower critical layer at 77 = 0, 
(where fj = A a Y — ac) and p 2 is an unknown function of £ and X. 


We now consider the Regions IV ± near to the upper critical layer situated at Y = Yu 
and write Y = Y\j + KY*. 


From the solutions in region II the flow quantities are expanded as 


u = /lAjy* + h 2 A*y* 2 + 6(uJ + hul + . . .), (3.8a) 

V = S{h 2 v{ + h\* 2 +...), (3.86) 

P — Pb + 6(/ip* + h 2 p 2 + ....), (3.8c) 

P = R* 0 + hRlY m + h 2 R;Y * 2 + 6(pJ + hp* + . . .), (3.8d) 
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w = s; + hSIY* + . . . + *K + hw 2 + • • •)> (3.8c) 

and substituting into the governing equations (2.1) yields at leading order 

vi = -£j= - A.,W - ac), (3.9a) 

^1*^0 

Pi = ( P i ex P(*7C) + c - c )- ( 3 - 9t ) 

At next order 

, 7 ~ \ . i r 2A 2« c i>I< , \ 

v 2 - -(^o f + acA 0x ) + A *^. \ a .2 + x*R^ P2< Plx ) 

+ H" 0 * 1 + * }) ' (31 ° a) 

P;=P 2 *(C,X), (3.106) 


where rfr = A*Y* - ac and the <f>f terms are from the continuation across the upper critical 
layer at fjx = 0, P* is an unknown function of X and p| is an unknown function of ( and 
X. 


The boundary layer flow must also match with the core flow outside of the boundary 
layer region hence we consider a far field region (region III Figure 2) where Y = Y /h and 
expand in the following manner 


u = Uoo + 6(hu i + h 2 u 2 + ....), 

(3.11a) 

v — 6(hvi + h 2 v 2 + ••••)» 

(3.116) 

p = Poo + S(hpi + ), 

(3.11c) 

P — Pb + 6(hpi + h 2 p 2 + ■••)• 

(3. lid) 

Substitution into the governing equations yields at leading order the Prandtl-Glauert equa- 
tion for pi with solution 

Pi = (Puexp(-7fiy)exp(t70 + c.c), 

(3.12) 


where 


and 



(3.13) 


^ = vi 0 exp(-7fty)exp(t7C), 


(3.14a) 
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t»10 = 


iP\ \ fl 

Poo Uoo 


(3.146) 


Real quajitities axe assumed and P\\ is an unknown function of X. At the next order 
the e* 7 ^ component of p 2 , say p 2 i, satisfies 


him + ^ 2 P2i (C = — 2ac£fooj^-pi a - Sl 2 pi X( - i-y^Pix- 


The solution is given by 


where 


P 21 = -P2iexp(-Tfiy)exp(i7C), 

ac(l-ft 2 ) 7 Piif , 


P 21 = P . 21 — 


siu 0 




(3.15a) 


(3.156) 


(3.15c) 


and P 21 (X) is a constant of integration. We note that higher harmonics are also present 
in p 2 . The e* 7 ^ component of the second order normal velocity disturbance v 2 can then be 
determined such that 

v 2 = v 2 oexp( — yilY)exp(ijC), (3.16a) 


where 


• rr . acjP n ' o5 
IP ooT^oo^20 — jj d” T**-* 21* 


(3.166) 


4 Derivation of the Amplitude Equation 

The am plitude equation governing the evolution of the mode is obtained from match- 
ing the expressions obtained for the pressure function and the normal velocity across the 
various regions of the flow. 

It is convenient to denote the e* 7 ^ components of P 2 ,p 2 , A^ , Aq , Ao <f> 2 by P 2 i > P 21 

respectively. The matc h i n g of pressure components across regions 11+ 

and III (equations (3.2c), (3.3a), (3.12) and (3.15)) yields at leading order 

Pu = Pi 1, ( 4 -l) 

and at next order matching gives 

P 21 = P 2 i - A^h, (4-2) 

where 

h = r ( p B 0% - !>~Vl)dY. (4.3) 

Jo 

In addition from the matching of the leading order normal velocity components in regions 
II * and III we obtain 

ftPn = AUl'rpoo. (4.4) 
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Matching the pressures between J"*’ and II shows that 

Pu = Pi 

and hence combining (3.6), (4.4) gives the dispersion relation 

IPooUlo = Aift-Roac. (4.5) 

Next matching the normal velocity components in regions II + and III at higher order 
gives using (3.16b), 

-At + Dii = . 2 a — P2u ( 4 - 6 ) 


where 


-KU.+lhA-rtf—n*, 

27 Yo iac ,Ai \ . 

Dl = |* 7Cr ~ J 2 _ p^Uoo + Tiu^\ ( X ) iiac u 


and 1%° is a finite quantity given by 

I°° - /°°I - - — \dY. (4.7) 

h " k W| PooUV 

Matching the pressure terms (equations (3.2c), (3.3a), (3.5c) and (3.76)) across regions 
1+ and II— gives 

Pn = Piu ( 4 - 8a ) 

and at next order 


Mat ching the normal velocity components across these two regions (equations (3.2a), 
(3.36), (3.56) and (3.7a)) yields 


—i'yacA^y-^- + = — 


where 




1 ,-RJ 2AJ . 1 . , Yu , 1 i n \Y\( R * I 2A *V1 

R 0 \lY 0 + (ji; + AJ ^JRJAJ 2 Fo-Fu iToAf 0 R 0 Aj 


j(°) _ fr* L_ 

2 ~/*W| RIKHY- 


IY 0 

1 


y u )2 RoXlY 2 


5 +^)+lr^(£+^)]^ 


, i r— 1. 4- 4 (— 4 : 

+ R;A; ! (y - -R; xi 1 R*xiy v r> x l 
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Further on applying 
velocity component at the 


the boundary condition at the wall (namely that the normal 
wall is zero) (3.7a) implies that 


-i^-2acix-^ + ,T “ 

it 0 Ai 


+ +ln 


.. , i~fR\a 2 c 2 - . 

ac\)A-\ — A = 0. (4.8d) 


J?oAi 


We now consider matching the normal velocity components (equations (3.2a), (3.36), 
(3.9a) and (3.10a)) across regions II and IV. This shows that 


A+ A- iljjLl f 2 ^2 I 

An Au xfR; l Aj + J^ 21 ^ !1 ' 


(4.9) 


Then on e limi nation of A 11 and between (4.6), (4.8c), (4.8d) and (4.9) leads to 
the amplitude equation 


DzA — 2 occAx — 




(4.10) 


where Dz is a purely imaginary constant defined by 


D 3 


'tUooPco „ , i~lR\{ac ) 2 , i~r(ac ) 2 

H^n Dl - nr + R „A, + A, 



jULpc* d 
R o\ 2 n 2 ' 


Equation (4.10) shows that the amplitude of the cross-flow vortex is dependent on the 
jumps ($ 1 ! — <£~ : ) and (^ — <f>2i) across the lower and upper critical layers. We consider 
next the details of the critical layers. 


5 Critical layer Analysis 

By increasing the size of 8 it is found that when S is of order 0(h 3 ) nonlinearity 
significantly alters the linear growth rate. This is the same as in Gajjar (1994) where 
some further justification and a balancing argument may be found. It is interesting to 
note that in this long wavelength limit nonlinearity gives rise to an unsteady nonlinear 
partial differential equation system governing the evolution of the disturbances, in contrast 
to shorter waves where nonlinearity leads to a single Hickernell (1984) type of integro- 
differential equation. 

We first consider the lower critical layer region (of thickness 0(h 2 )) and write Y = 
h + h 2 fj. The flow quantities axe expanded in the following manner 

u = hac + h 2 (x 1 rj + ^r~) + h 3 U* + h*U; (5.1a) 
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(5.16) 

(5.1c) 

(5. Id) 



The upper critical layer, at Yj/, is considered next and we write Y 
and the flow quantities there expand as 

ky + /iff +/i 2 77i 
A 1 

a = hac + + k’or + . 

(5.5a) 

v = h 5 v* + tfv; + h 7 V : * + ...., 

(5.56) 

p — pb ■+■ h* P* + h? P 2 d" • • • > 

(5.5c) 

P = R 0 + A"" + *i + • • • • 

A i 

(5.5d) 

Substitution into the compressible Navier- Stokes equations yields 


L i($i)=0, 

(5.6) 

Ll{qU) — p* 2 5 

(5.7) 
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where L\ is the operator given by 


and 



<y 4- 4- 


XI a 2 c 2 'i 

X * 2 J 




+ ac 


d 

dX 



(5.8a) 


Vo* = 


P* 


x\K 


QU 


= P£h- 


(5.8i>) 


The solutions in the lower critical layer region must be matched with the solutions 
in regions J ± (of thickness 0(h)) just outside of the critical layer so as to determine the 
behaviour of Ql and Hi as rj -» ±oo. Matching using (3.5) and (3.7) yields 


q L ~ 2X 2 fj 4- + ^•}- Re (^ ex P(*'70) + ••• 


as 


Tj -¥ ±00, 


(5.9) 


and 

Hi ~ Rirj+ — r^-r--Re(Aexp(t7C)) + ••• as — ► ±oo. (5.10) 

XiV 

Similarly from matc hin g the upper critical layer with regions IV equations (3.9) 
and (3.10) yields, 


and 


as 


1 U ~ 2 A ^1 + 2 ” c a A .^° { } M* «p (nO) + • • • > 

$1 ~ + 2 -Re(X exp (i T C)) + • • • , 


(5.11) 


(5.12) 


77l “4 ±00. 

The solutions of (5.2), (5.3), (5.6) and (5.7) together with the amplitude equation, (4.10), 
govern the evolution of the vortex. We may rewrite the jump conditions in the form 



qL exp (-i'tC) drj dC, 


2irac 

7 



(5.13 a) 


and 


Jj y’“ we rp(-nO<% <«C = (^ + - «■)• ( 5 - 134 ) 


with f* denoting a finite part integral. This enables us to rewrite the amplitude equation 
as 
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q L exp (-ij()drjd( = 


- t7 2 ac f -* f*°° 

D 3 A - 2 acAjc — - . / / 

ZTTAi J 0 J-oo 

/ 7 / ?uexp(-z7C)<i»?i<*C> 

27TAJ Jo J-oo 

where the .D 3 term gives a correction to the wavenumber. 


(5.14) 


A more standard form of the equations is obtained by re-nor m a li sing using tbe fol- 
lowing scalings: 

A 2 aV , 

Aj 77 + — r 0 — Aid^Zi, 


AJ 771 + 


Ai 

A;<»V 

K 


+ 0 — AjdJ^i, 


Hi = Rid*n, $1 = 


jr=-^-l-x 0 , tC = C + Co, 

Ai7“* 

A = dJJAexp (-»Co), 
qL = ^Ql, qu = d$Qu, 

(where Xo and Co are chosen so as to match with upstream conditions). The equations 
(5.3), (5.2), (5.7) and (5.6) respectively are rescaled as 


^ + Z a AL - MUexp (O)^y- = JRe(iAcxp(iC))^, 

an .an ...... an 

= JiJt*Aexp(0)||-. 

and the amplitude equation (5.14) is transformed to 

.2ir *00 


(5.15a) 

(5.15fc) 

(5.15c) 

(5.15d) 


Ai = - i h. r [°° Q L txp(-iC)dCdZ- 

* Jo J-00 


^ r r Quexp (—i(*)d(*dZi + D 3 A 

^ Jq J —00 


(5.15e) 
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where 


J = 


Ex 

Ho 




R’ 

H’ 


The constants d * , A,£i, E 2 are given by 

d{ = 


a 2 c 2 


2 \\ ’ 

Ai^ 2 


* = lac ’ 

2 2O!cAli^0^Q 


4r- 


a; 2 *; 


Dz = 2Ai'y<^23j 

(1- J)^«2A 2 <Q, 

(1 - h)dl = 2A^ 4 *, 


(5.16a, 6) 


,2A; h; fioAf 

a; + RVRIK 2 

Equations (5.15) must be solved subject to the following boundary conditions 




Ql ~ (1 — J)Z + 

n~ z + 

Qi/ ~ (1 — J\)Z\ + 
$ ~ Z\ + 


Re(Aexp (i£*)) 


fle(Aexp(t'C*)) 

Z 

He(Aexp (tC*)) 

Ti 

Re(i4exp(zC*)) 

Zi 


Z — ► ±oo, 

.Z — ► ±oo, 

, ^1 ~ ► ±00, 

Z\ —* ±oo, 


(5.17a) 

(5.176) 

(5.17c) 

(5.l7d) 


The norm aliz ed linear growth rate is K r (defined below) so that Re(A^/A)-^K r as X-*-oo. 

We note that the equations are now of a similar form to the evolution equations 
obtained by Goldstein and Wundrow (1990), and Gajjar (1994). The equations (5.15) 
are a generalization of those of Gajjar (1994) for flows with multiple critical layers. In 
our analysis however viscous effects have been ignored (such effects come into play when 
h = 0(H) -1 / 14 ), see Gajjar (1994). The equations are then modified to 


CQl ~ 1cQ L zz = J Re ( iA «P(*C*)) n # + 7«n**, (5.18a) 

ni-7<*n^ = o, (5.186) 

CuQu -liQuitix - J i Re ( iAex P( i <*))*2 1 + (5.18c) 

j Cu$ — 7d$z 1 £ 1 = (5.18d) 
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where the operators C^Cu a-nd the constants in (5.18) are given by 


c ^H + ^- Re(iA ' Min) £' 
Lu e +Z 'W- - 


SoniS^R-’h- 7 

7c " , M; 3 7 

Prix'(S 0 ) . 

= 7 - J( " } ’ 


-w = Tc/-P r > 

7e n(So)R;K< 


a=ti!Pr, 7.* = 7jJi - D- 


(5.19) 


Here /x(T B ) is the viscosity law (5 0 is defined in (2.4)), and because of the assumptions 

made in section 2, 7 C> 7* take opposite signs. 

The equations (5.18), with boundary conditions (5.17) govern the nonlinear evolution 
of the cross-flow vortex. Note that for a given basic flow the constants J, J\ are fixed and 
the only free parameters are 7 C and the wave-angle 6. For the inviscid case j c is zero. 
Linear Critical Layer 

The lin ear theory may be recovered from the above equations by first writing 


Ql =( 1 - J)Z + Re{Q L e iC ), 

H = Z + Re{He iC ), 
and then linearising the equations to give 


Qu = (1 — Ji)Zi + ReiQu^ )> 
$ = Z x + He($e iC * ), 


QlX +iZQL — IcQlZZ — + 7 e^zz' 

A Qux "b iZlQu — JcQuZiZi = ^ + le^Z 1 Zi’ 

The solution of these equations may be obtained by utilising the Fourier Transform tech- 
nique of Hickernell (1984). This gives 

fl f = iA(X + s)e‘ ,i * S/3 H(-ssga(jd))sgn('rd), 


= iA( X + As)e 7 ^ S/3 5(-ssgn(7a))sgn(7a) } 


Ql = ^ 


/ 3 iA(i" + s){l- 


37c /, ,7a-7«.s 

(1 — exp( — 5 


{id ~ 7c) 


))}b(-s sgn( 7 c )) sgn( 7 c ), 

(5.20a - d) 
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= e T ** 3 / 3 iA(X + As){l - -«P( ~ 3 7 -^ 3 ))}g(-^ s gn(7c))sgn(7e ) 1 


where the Fourier Transforms are defined by 


U^(s) = J U exp(-isZ) dZ, U* ex P(j sZ ) ds > 

with Z being Z, or Z\ as appropriate. Also in (5.20) fT(x) denotes the Heaviside unit 
step function. The above solutions agree with those given in Leib (1991) for s imil ar linear 
equations. The solutions (5.20) show that 


i a2tt a + oo ^ i /» 27 T a+oo 

- / Qlc~ 1<: ' dC dZ = i7TSgn(7c), — / / Que~ x ^ d(* dZ x = Z 7 rsgn( 7 *). 

T* Jo J—oo ™ J° J—oo 

6. Additional Comments and Discussion 

We may use the results in the previous section, or alternatively (5.13), (5.14) with 
the appropriate linear jumps, to obtain an expression for the linear growth rate of the 
cross-flow vortex. The scaled linear growth rate K r is thus given by 


.R o Ai 7 C 7 tcos 0 r 1 . 2 A 2 Ri 


[>k {2 -£ + - l )sgn(Ai) + + ! )sgnW) ] • (M) 


The expression (6.1) may also be expressed as 


RqXiTitccos6 / D (pbDUb) \ , nr - r J 
Kr_ 2 ^i( p 2 B D(U B y ) Sgn ^ Ub ^\' 


( 6 . 2 ) 


where the su mm ation is taken over all critical levels and the expressions inside the sum- 
mation axe evaluated at the critical levels, and D = 

As in Gajjar (1994) the unsealed wavenumber, frequency and growth rate terms are: 

wavenumber R* 7 * = R 7 {~fh + ...), frequency R 7 = R 2 (jh 2 c cos 8 + — ), 

growthrate R*i c* = R 7 (K r h 2 +...). 

Using (3. 13), (4.5) the leading term in the expansion for the wavenumber is given by 


-£ = 1 - v&il)*- 

^ OO 


The result (6.2) shows that for neutral waves we require 
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which is a generalization of the result that D(p B DU B ) = 0 for flows with a single critical 
level. 

Unlike the planar case the result (6.1) shows a more complex dependence of the 
growth rate on the wave-angle 6 since Aj, AJ,A 2 ,AJ! depend on $ through (2.5), (2.6). For 
the neutral case the condition that U B = 0 together with (6.3) fixes the position of the 

critical level as well as 9. « 

The second important observation ste mm ing from (6.1) is as follows. For a two- 
dimensional insulated, flat plate boundary layer with zero pressure gradient, the growth 
rate k t is zero, Gajjar (1994). For a fully three-dimensional boundary layer this is no 
longer true because of the presence of the last two terms in (6.1). This shows that for the 
stability problem with cross-flow effects present, there is no real distinction between the 
insulated and non-insulated cases as the growth rates are of the same order in both cases. 

The analysis used to derive (6.1) was based mostly on an inviscid analysis with viscous 
effects confined to the critical layers. The Stokes wall layer also contributes to the growth 
rate and as Appendix A shows this contribution can be written as 


R * /c* = R 2 K r h 2 + • • • + R* h 


i 7*1 
2\/2|c| cos0 


( MSo) \* 

\ilo|c|7cos0/ 


+ 


(6.4) 


where the last term in (6.4) arises from the viscous wall layer. For 9 not close to f and 
h » R~ the Stokes layer contribution in (6.4) is much smaller than the inviscid con- 
tribution. As h decreases however the two terms become comparable when h 0(i? 12 ). 
This implies an effective wavenumber of O(RTi) and growth rates of O(-Rs). This is pre- 
cisely the scaling for the upper-branch of a two-dimensional accelerating boundary layer 
in incompressible flow considered by Smith & Bodonyi (1982), and extended to the three- 
dimensional boundary layer case by Bassom & Gajjar (1988). In this limit it is noted 
from (5.19) that 7 C ,7* are large parameters and as Appendix B shows, for | 7 C | » 1 the 
unsteady analysis remains largely intact but with the critical layers beco min g strongly 
viscous critical layers of the Haberman (1972) type. Additional diffusion layers at the 
outer edges of the critical layers are necessary to reduce the large mean vorticity jump 
produced by the critical layers back to zero outside in the main part of the boundary layer. 
When h = 0{R~&), from (5.19) it can be seen that 7c = 0(R&). From Appendix B the 
thickness of the critical layer is then of 0(R~2h 2 yi ) = 0(R as) which is the same as in 
Smith & Bodonyi (1982). The diffusion layers at the outer edge of the critical layer have 

thickness 0(i? _2 /i 2 7c ) = 0(H~s6). 
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Hence in the formal limit when h becomes of 0(R « ) the linear growth rate is given 


by 


KU. = Ri + l^gntA.) + + §)sgn(A;)l 

\ Z LAj/to Aj -fto Aj xto Aj JTLq J 


+ 


7*1 

2\/2|c| cos# 


/ Mgo) 

\ilo|c|7cos^/ y 


(6.5) 

The condition for neutral waves is obtained by setting the expression given in (6.5) to zero 
and it can be seen that a special case of this result agrees with that obtained by Bassom &: 
Gajjar (1988) for incompressible flow. (The incompressible results are retrieved by setting 
PB to be constant). 

When h = 0(f?~“) results for the unsteady nonlinear evolution of the cross-flow 
vortex may be s imil arly deduced. The appropriate expression for the leading term of the 
growth rate is then, 


R 


*(- 


7Xld *-[EjTn(i<f> L ) + E 2 Im{i<f>u)] + 


ac 


7*i 

2\/2|c| cos# 


(4!Sb)> ■ 


where the nonlinear jumps <f>L , 4>u a t the lower and upper critical layers respectively, are 
obtained from 

1 -2* »*oo «i p2tc *+oo 

A<f>L = - / e~ iC Q L dZdC, A<f>u = — / / e-^QudZdC, (6.7) 

JO •/ — oo Jo J — oo 


and Ql,Qu satisfy the steady versions of (5.19), (see also appendix B). Some results for 
<f>L for the single critical layer version of (5.19) ancLwith the Chapman viscosity law are 
presented in Gajjar & Cole (1989). For nonlinear neutral waves the expression in (6.6) is 
set to zero and this again extends the Bassom & Gajjar (1988) results to the compressible 
case. 

The full unsteady nonlinear problem governed by (5.17), (5.19) involving the solution 
of four coupled unsteady nonlinear critical layer equations poses a formidable numerical 
task and is not pursued here. The problem is being tackled using the methods, described 
in Gajjar (1994), for the single critical level case. 
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Summary 

In conclusion we have obtained the nonlinear equations governing the evolution of 
a long wavelength cross-flow vortex in a three-dimensional compressible boundary layer. 
These equations remain to be solved but the solution properties axe of some interest espe- 
cially in relation to the question of whether the cross-flow vortices attain an equilibrium 
amplitude or not. The one major s im plifying assumption used in the present work is that 
the cross- Squire variations are small and negligible. The relaxation of this assumption 
requires further study but it is expected that as the z variations gradually increase the 
first major non-trivial influence will be felt in the critical layers. The question of shorter 
waves is addressed in Gajjar (1995) for the stationary cross-flow vortex where a novel 
integro- differential amplitude equation for the vortex has been obtained, and in Gajjar & 
Arebi (1995) for the non- stationary case. 

Although it has not been done here, it would be interesting to compare the predic- 
tions for the linear growth rate as given by (6. 1-6. 5) with full numerical solution of the 
compressible linear stability equations. Comparisons for the single critical layer case are 
shown in Gajjar (1994) and good agreement was found particularly for the insulated wall 
case. A special case of (6. 3), (6. 4) includes that for a two-dimensional boundary layer and 
it would be interesting to see how the viscous correction to the growth rate compares with 
the equivalent £ Orr-Sommerfeld 5 type or linear stability calculation for the first mode. The 
results (6.3),(6.4) obtained here are particularly useful as they provide an alternative route 
for investigating the effects of different parameters such as Mach number, wall heating, etc 
on the stability properties. U nlik e the numerical solution of the linear stability equations, 
the analytic results given here may be used in a much more systematic and easier manner 
for investigating these different effects. 

Finally we have shown how the unsteady critical layer analysis may be linked to the 
properties of the steady Haberman type of critical layer. In the limit that the Haberman 
par am eter 35 used in this work) becomes large, the major modification required in 
the unsteady critical layer structure is the inclusion of appropriate diffusion layers at the 
outer edges of the critical layers. The analysis in Appendix B apply equally to other flows 
including the accelerating boundary layer and channel flow as considered for example by 
Gajjar & Smith (1985). An application of such ideas to channel flow with compliant 
boundaries is presented in Gajjar &; Sibanda (1995). 
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Appendix A 


In this appendix we derive the correction to the inviscid growth rate arising from the 
Stokes wall layer. For the sake of clarity only the important details are given. 

We first note that the thickness of the Stokes layer is given by 

y = R-*h~ 1 y , or Y = R~*h~ 2 y, (Al) 


from a balance of the viscous and inertial terms in the Navier-Stokes equations. For the 
layer to be thin we require that ft » 1?~«. Expansions for the flow quantities are readily 
written down using (3.4), (3.5) and in particular we have 

u = R 4 ft. 1 Aiy + Sui + . . . , 
v = SR 4 + . . . , 

w — 6 2 R~* h _1 y + 6R~*h~ 2 wi + . .. , (A2) 

P = PB + Shpi + . . . , 

p = JZo + SR 4 h 2 pi + 


Substitution into the Navier-Stokes equations yields the usual equations for the Stokes 
wall layer. These are easily solved and the finite part of the normal velocity component 
as y — » oo, say V\ is calculated to be 

_ Pic _ ^i^QC 
Vl f macRo m ’ 



where 


V M 5 o) / 

This shows that the expansion for the normal velocity in zone I± is of the form 

V = ^(ft 2 U! + fl*V 2 + . . . + R + •••)> 


(A4) 


where we require 

v ti {Y = 0) = v lf . (AS) 

To find the effect of the displacement (A3) on the growth rate it is necessary to 
introduce an additional multiple- scale coordinate X 2 defined by 


d 

dx 




d 

dX 2 ’ 
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I - + ^ + “ 4 > + 

where u is a constant to be found. We now regard Aq as a function of X 2 <dso so that 
Aq = (A(X,X 2 )exp(i() + c.c). 

The expansions in the various regions are also modified so that, for example, in the 
Y = 0(1) region we have 

u = Ub + + hu 2 + . . . + R 4 h?u s i + ...), 

v = <5(/itq + h* v 2 + . . . + R 4 h ^ v s i + . . ■)> 

and in region I± we have 

u = /iAjF + . . . + f(Si + . . - + 22“ + ...), 
v = 5/i 2 (£i + . . . + 22“ 4 Vji + ...,) 

p =: PB + 8(hpl + . . . R 4 h 1 p3t+--*) 

For the solutions to match we require v = — 2. Substitution of (A6),(A7), and similar 
expansions in the other regions, into the Navier-Stokes equations leads to a sequence of 
linear equations for the ‘sz’ quantities which are easily solved. For example 

v,i = p~T < PsiC + P1X2) - T~ UlX l - ~ ( A8 ) 

KqAi Ai 

and so on. Solving for each of the i sV components in turn and matching then leads to 

- e^Mk 
Psi ~ n 

Hence setting Y = 0 in (A8) and using (A4),(A5),(4.4),(4.5) shows that 

***> = ~2 

This gives the correction to the growth rate as stated in (6.4). 


(A6) 


(A7) 
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Appendix B 

In this appendix we consider the solution of the system (5.18) with boundary condi- 
tions (5.16) in the limi t a s |-y c | becomes large. Without loss of generality we will take j c to 
be positive and write 7 e = eiT’cj where ej is a constant. It is sufficient to consider just the 
lower critical layer as the analysis for the upper critical layer follows analogously. Hence 
the governing equations (with some of the suffices dropped) are 


d Q; + Z |p- - fle(iAexp(»0)|| - 7 c ^7 = JRe{iAexp(iC))j£ + e i 7 c|^, (Bio) 


dX 


du.an .... a n lc a 2 u 

+ Z— - Re(iAex p(t< ))^ - — — = 0. 


az PrdZ 2 


(B16) 


with boundary conditions 


q ~ (i — j)z + + ... Z ->± oo. 


n~z + JU(At ^ (iC>) +... z~± oo. 


(B 2 o) 

(B26) 


For 7 C >> 1 there are two distinct regions which occur, see Figure 4. A balance of 
ZQf term with the viscous term shows that Z = 0(j c s ). This region, which we will denote 
by CLl, comprises the main part of the critical layer in which the solution properties are 
similar to those for a Haber man type of steady (or equilibrium) critical layer. A balance of 
the Qx term with the viscous term shows that Z = 0 ( 7 ,? ), and we will denote this region 
by the label Dl. One of the important properties of the Haberman critical layer in region 
CLl is that it generates a mean vorticity jump, see for instance Stewartson (1981). The 
main function of the region Dl here, or diffusion layer to use a more fa mili ar terminology, 

is to restore the mean velocity jump to zero outside the critical layer as required by (B 2 ). 

12 - 

Thus in region CLl we set Z = yi Z u A = ~fi A and we have the expansions 


Q = Qo + • • • > H = Ho + — . 


Substitution into (Bl),(B 2 ) then shows that at leading order Qo ; Go satisfy 


* 7 $ - ^M«p(.C))|| - § = 


(B3a) 

(B36) 
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These are just the compressible version of the Haberman equations, see Gajjar & Cole 
(1989), and they have the property that 

0. ~ (1 - J)Z, + B±(X) + + . „ (B4fl) 

no~Z l + A ± (X) + ^ta^ + ... Zi - ± 00 . (B4t) 

The quantity (#+ - H~) is the mean vorticity jump produced by the critical layer. Note 
that in (B 3 ),(B 4 ), A,H ± ,A ± are all dependent on X, but since X appears only as a 
parameter these equations are solved exactly as in the steady case. 

It is convenient to treat the diffusion layers above and below the critical layers si- 
multaneously. Hence in region Dl^, see Figure 3, we set Z — 7 c V an d (B4) imply the 
expansions a s 

Q = 7 c?o + 9 i + 7<T*92 + 7c 6 93 + 7c 6 94 + ---> (B5ab) 

n = 7 / n 0 + Hi + 7 C 6 n 2 + 7 c 6 n 3 + 7 C 6 n 4 + — 

Substitution into (Bl) then yields a sequence of equations given by 

77 II 0 C = 0, 79 o{* '= 0 j (B 606 ) 

77 H 1C . = 0, Tjqif = 0, (B7ab) 

7 ?n 2< . - Re{iAe iC )Ho„ = 0, 77 q 2f - - Re(iAe iC )q 0v = JRe{iAe' c )H 0 „ (B8ab) 

t/Hs^* + n oJ £ - Re(iAe l( )Hi, — -p^^-Ovv = °> (B9ab) 

VQ3C + %x ~ Re{}Ae l( ' )qi v ~ 9otj? ? = JRe(iAe t( * )Hi, + eillo^r?) 

tfUf + - i?e(iAe lC )H 2n - = °» (BlOab) 

WQiC + ?! x ~ Re(iAe t< ' )q 2 17 — 9ivn = ^ Rt{iAe^ )H 2n + eilli^. 

These equations may be solved in series. First (B6),(B9) show that 

Ho =V, 9o = (1 - J)*l- ( B11 ) 

The solution of (B7) gives 

Hi = HimC-^j 7 ?)? 9i = 9i m(X,t]). (B12) 

The solution of (B 8 ) yields 

n 2 - R < A ' iC .i +n 2M (x,n), q2 = R ' (Aet< ) (B 13 ) 

T) V 
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If we substitute for Ho, qo, Hi, gj , II 2 , 92 into (BlO) and use (B11-B13), then after averaging 
(BIO) over a period in gives 


dlli M _ d 2 IIiM _ dqiM _ d 2 <7iM _ d 2 HiM 
dX dv 2 ~ ’ dX dv 2 ~ 61 dv 2 

To match with region CLl we require that 

n 1M (X, ± 0 ) = A ± (X), q 1M (X, ±0) = #±(1), 


(Bl4a, 6) 


(Bl4c, d ). 


The problem (Bl4a,c) is identical to that solved by Smith et a/.(1993) in their investiga- 
tion of vortex wave interactions. The solution of (B14) may be obtained using Fourier 
transforms and it cam be shown that the solutions which decay to zero at 77 = ±00 are 
given by 






1: 


(X - s)f 


A 


ds. 


If Pr = 1 then 


db 

9im 


£ 


JsL 

,2a/7T 


B*(s) 




ds 


- ex 





ds , 


and if Pr ^ 1 then 


qiM = 


I v 


2Vn 


r.[( 


H ± (s) + 


eiA^aJN exp ( 4(l-s) ) 


-1 )) 


(Pt-I)J (X-s)f 

nfPr 


e 1 A ± (s)\/Pr ex P( 

(Pr-1) (j£-s)§ 


ds. 


The ± refers to the solutions in the regions above and below the regions CLl respec- 
tively. 
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Figure 1.— U B > W B> are as in the text. Y u is the location of the 

upper critical layer, (a) A sketch showing the coordinate system 
with the Squire coordinates £ and z. (b) A sketch of typical 
velocity profiles for a three-dimensional boundary layer. 
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Figure 2. — A schematic diagram showing the various regions of the flow including the critical layers. 
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Figure 3.— A sketch of the diffusion layers DL* and the inner critical layer region CL1. 
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